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Ebert et al. [Phys. Rev. Lett. 77, 3827 (1996)] have fractured icosahedral Al-Mn-Pd single crystals 
in ultrahigh vacuum and have investigated the cleavage planes in-situ by scanning tunneling mi- 
croscopy (STM). Globular patterns in the STM-images were interpreted as clusters of atoms. These 
are significant structural units of quasicrystals. The experiments of Ebert et al. imply that they are 
also stable physical entities, a property controversially discussed currently. For a clarification we 
performed the first large scale fracture simulations on three-dimensional complex binary systems. 
We studied the propagation of mode I cracks in an icosahedral model quasicrystal by molecular 
dynamics techniques at low temperature. In particular we examined how the shape of the cleavage 
plane is influenced by the clusters inherent in the model and how it depends on the plane structure. 
Brittle fracture with no indication of dislocation activity is observed. The crack surfaces are rough 
on the scale of the clusters, but exhibit constant average heights for orientations perpendicular to 
high symmetry axes. From detailed analyses of the fractured samples we conclude that both, the 
plane structure and the clusters, strongly influence dynamic fracture in quasicrystals and that the 
clusters therefore have to be regarded as physical entities. 
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I. INTRODUCTION 

Quasicrystals are intermetallic compounds with long- 
range quasi-periodic translational order. They possess 
well-defined atomic planes and hence diffract electromag- 
netic and matter waves into sharp Bragg spots. But t hey 
also display atomic clusters as basic building blockalEI^ 
whose arrangement in space is compatible with the pla- 
nar structure. These clusters consist for example of sev- 
eral shells of icosahedral symmetry (Bergman-, Mackay-, 
pseudo-Mackay-clusters). Or they form polytopes, e.g. 
decagonal prisms, which like the unit cells of periodic 
crystals fill space, although with large over laps ("quasi- 
unit-cell picture" )0. Janot and otherd^*^ have postu- 
lated that a self-similar hierarchical assembly of the clus- 
ters is responsible for the stability of quasicrystals and 
for many physical properties, like the low electric conduc- 
tivity. However, it is a controversial and persistent dis- 
cussion, whether the clusters are merely structural units 
or whether they represent physical entities. The discus- 
sion was fueled by an experiment of Ebert et al.^^, where 
icosahedral Al-Mn-Pd was fractured under ultrahigh vac- 
uum conditions at room temperature. Scanning tunnel- 
ing microscopy images of the cleavage planes revealed 
elements of 0.6 to 1 nm in diameteif^. The authors ar- 
gue that these are the clusters which were circumvented 
by the crack and hence form highly stable aggregates of 
matter. Others point out that flat terraces evolve on 
fivefold surfaces of i-Al-Mn-Pd when annealed at high 



temperatured^Hnmnil, As this requires truncated rows of 
Bergman and Mackay clusters it is stated-'^ that these 
therefore could not represent firm entities. 

In the present article we report on molecular dynamics 
simulations of crack propagation in a three-dimensional 
icosahedral model quasicrystal at low temperature. Seed 
cracks are inserted along different planes and therein 
along different directions. The fracture planes are care- 
fully analysed to answer the role of clusters in dynamic 
fracture. In Sec.|ll]we provide some requirements on the 
theoretical description of fracture. In Sec. |III| the model 
quasicrystal, the molecular dynamics technique, and the 
methods to visualise the results of the simulations are 
outlined. Subsequently, in Sec.|IV]the simulation results 
are presented and then discussed in Sec. |V] 



II. FRACTURE 

The stress concentration and the strength of the load- 
ing at a crack tip are determined by the macroscopic 
geometry and dimensions of a sample. In linear elastic 
continuum mechanics, a sharp mode I (opening mode) 
crack is characterised by a singular stress field and a cor- 
responding displacement field, which both are propor- 
tional to the stress intensity factor K. This factor is 
proportional to the applied external load and contains 
the geometry of the sample. A simple energy based con- 
dition for crack propagation is the Griffith criteriorfSl. 
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It states that a crack is in equilibrium when the change 
in mechanical energy per unit area of crack advance - 
the energy release rate G - equals the change in surface 
energy of the two fracture surfaces, 27. In continuum 
mechanics the energy release rate is proportional to the 
square of the stress intensity factor for a given mode. A 
crack then should start moving when the stress intensity 
factor exceeds the critical Griffith value. 

A continuum mechanical description of fracture, how- 
ever, has a few drawbacks. First, the requirements for 
linear elasticity are no longer valid near the crack tip 
where atomic bonds clearly become non-linear and even- 
tually break. Second, a continuum theory neglects the 
discrete nature of the lattice. Thus, it is fully ignored 
that fracture of materials is ultimately caused by bond 
breaking processes on the atomic scale. 

A way to understand the processes is to perform nu- 
merical experiments, since experimental information on 
this length scale is difficult to obtain. Molecular dynam- 
ics studies have provided useful insight into crack prop- 
agation in pure metals and simple intermetallic alloys, 
whereas in complex metallic alloys the mechanisms are 
not yet so clear. Atomistic studies show for example 
that cracks remain stable in a region around the critical 
stress intensity factor due to the discrete nature of the 
lattice. This effect is called lattice trapping'^. A further 
consequence of the discrete lattice is that the fracture 
behaviour in one and the same plane can depend on the 
crack propagation directio^^. Such observations cannot 
be explained by a simple continuum mechanical descrip- 
tion. 



III. MODEL AND METHOD 

A. Icosahedral binary model 

In the numerical experiments we use a three- 
dimensional model quasicrystal which has been proposed 
by Henley and Elsei'^ as a structure model for the icosa- 
hedral phase of (Al,Zn)Mg. This is the simplest possi- 
ble model quasicrystal that is stabilised by pair poten- 
tials. Furthermore it allows Burgers circuit analysis and 
is a prototype of Bergman-type quasicrystals. As we 
do not distinguish between Al and Zn atoms, we term 
this decoration icosahedral binary model. It can be ob- 
tained by decorating the structure elements of the three- 
dimensional Penrose tiling, the oblate and the prolate 
rhombohedra (see Fig. [l] top). Al and Zn atoms (A 
atoms) are placed on the vertices and the midpoints of 
the edges of the rhombohedra. Two Mg atoms (B atoms) 
divide the long body diagonal of each prolate rhombohe- 
dron in ratios t:1:t, where r is the golden mean. Two 
prolate and two oblate rhombohedr a wit h a common ver- 
tex form a rhombic dodecahedroEp^*^. To obtain the 
icosahedral binary model, in these dodecahedra the atom 
at the common vertex is removed and the four neigh- 
bouring A atoms are transformed into B atoms. Finally 



these atoms are shifted to the common vertex to divide 
the edges of the corresponding rhombohedra in a ratio 
of 1:t. Fig. [T| (bottom left) shows the final decoration of 
these dodecahedra, in which the B atoms form hexagonal 
bipyramides. This modification increases the number of 
Bergman-type clusters (see Fig. [l] bottom right) inher- 
ent in the structure, leads to a higher stability with the 
potentials used, and takes better into account the experi- 
mentally observed stoichiometry of the quasicrystal. The 
Bergman-type clusters may also be interpreted as build- 
ing units of the quasicrystal and are the main feature of 
the structure apart from the plane structure. 



B. Molecular dynamics technique 

As is often done in fracture simulations (see e.g. Abra- 
hampfll) we use simple Lennard- Jones potentials to model 
the interactions. The following facts led to this choice: 

First, the very few potentials available for quasicrys- 
tals are unsuitable for fracture simulations: These pair 
potentials stabilise the quasicrystal only for a fixed vol- 
ume. With free surfaces atoms sometimes simply evap- 
orate (e.g. with potentials based on those of Hafner et 
al.*^. Very long-ranged potentials with Friedel oscilla- 
tions (e.g. those of Al-Lehyani et al.l^ frequently display 
such large cohesive energies that nearly no elastic defor- 
mation is possible and instead intrinsic cracks develop. 

Second, many simulations have proven that model po- 
tentials are helping to understand the elementary pro- 




FIG. 1: Tiles of the icosahedral binary model decorated with 
two types of atoms and the Bergman-type cluster. Top: pro- 
late rhombohedron (left) and oblate rhombohedron (right). 
Bottom: rhombic dodecahedron (left) and the 45 atoms build- 
ing the Bergman-type cluster (right) inherent in the model 
quasicrystal. Small grey and large black spheres denote A 
and B atoms respectively. 
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cesses in fracture (see e.g. AbraharrP^ and so are rea- 
sonable when qualitative mechanisms are the centre of 
interest. For quantitative results, a known limitation is 
the neglect of non-local and many-body interactions. 

Third, the Lennard-Jones potentials use cpl23keep the 
model stable even under strongest mechanical deforma- 
tions or irradiation (introduction of point defects) and 
have been used in our group in many simulations of dis- 
location motioii^ or even shock wave^^U^ 'pjjg structure 
is robust under a wide variation of the potential depths. 
Very similar potentials have shown to stabilise the icosa- 
hedral atomic structure in a simpler modeP^. it is also 
known since the early fifties that Lennard-Jones poten- 
tials favour icosahedral clusterd^, indicating that these 
potentials are useful for structures like icosahedral qua- 
sicrystals. 

The minima of the potentials for interactions between 
atoms of the same type are set to eqi whereas the mini- 
mum of the potential for the interactions between atoms 
of different kind is set to 2eo. The conclusions drawn 
from our simulations, however, remained essentially un- 
changed if all binding energies are set equal. The shortest 
distance between two A atoms is denoted tq. All masses 
are set to mo. The time is then measured in units of 
to = roy^mo/eo. 

The molecular dynamics simulations were carried out 
using the micr ocanonical ensemble with the program 
code IMEP^MI, Jt performs well on a large variety of 
hardware, including single and dual processor worksta- 
tions and massively parallel supercomputers. 

First, we searched for the potential cleavage planes. 
According to the Griffith criterion they should be planes 
of low surface energjPi'. To identify these surfaces we re- 
lax a specimen and split it into two parts. Subsequently, 
the two regions are shifted apart rigidly. The surface en- 
ergy is then calculated from the energy difference of the 
artificially cleaved and the undisturbed specimen. 

In contrast to simple periodic crystals, the atomistic 
structure of the planes and therefore also the surface en- 
ergy in quasicrystals varies strongly within the family 
of planes perpendicular to a fixed axis. In the present 
model, surfaces with low interface energy are located 
perpendicular to two- and fivefold directions at certain 
heights (Fig. |2|. Perpendicular to other directions the 
plane structure is less pronounced and the minimal sur- 
face energies are higher. 

Since we are interested in the morphology of fracture 
surfaces we apply a sample form that allows us to fol- 
low the dynamics of the running crack for a long time. 
For this purpose, a strip geometry is used to model crack 
propagation with constant energy release rat^^. The 
samples consist of about 4 to 5 million atoms, with di- 
mensions of approximately 450ro x ISOrg x TOrg. Periodic 
boundary conditiond^ are applied in the direction paral- 
lel to the crack front. For the other directions, all atoms 
in the outermost boundary layers of width 2.5ro are held 
fixed. An atomically sharp seed crack is inserted at a 
plane of lowest surface energy, from one side to about 
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FIG. 2: Surface energy of cleavage planes perpendicular to 
two- and fivefold axes. 

one quarter of the strip length. The system is uniaxially 
strained perpendicular to the crack plane up to the Grif- 
fith load and is relaxed to obtain the displacement field of 
a stable crack at zero temperature. Then a temperature 
of about 10~^ of the melting temperature is appliecP^to 
the configurations with and without the relaxed crack. 
From the resulting configurations we obtain an averaged 
displacement field for this temperature. The crack now is 
further loaded by linear scaling of this displacement field. 
The response of the system then is monitored by molec- 
ular dynamics techniques. The sound waves emitted by 
the propagating crack (see Fig. [s] and online movid^ 
are damped away outside of an elliptical stadiurrPH to 
prevent reflections. 

C. Visualisation 

To study crack propagation on an atomic level the se- 
lection and reduction of data is of crucial importance. 
Due to the large number of atoms required for the study 
of crack propagation in three-dimensional systems, it is 
not feasible to write out the positions of all atoms during 
the simulation, and even less to display all of them. To 
obtain a first insight into crack propagation only atoms 




FIG. 3: Kinetic energy density. Sound waves emitted by the 
propagating crack and the elliptical region without damping 
are clearly visible. See online movie^. 



FIG. 4: Snapshot of a simulation with some 4 milhon atoms. 
Only atoms with low coordination number are displayed. See 
online movi^^. 



FIG. 5; Bergman-type cluster cut by a flat surface. 



near the fracture surfaces are of interest. Whereas they 
can be visuahsed in periodic crystals by plotting only 
those atoms whose potential energy exceeds a certain 
threshold, this technique is not applicable for quasicrys- 
tals. Because of the largely varying environments the 
potential energy varies significantly from atom to atom, 
even for atoms of the same type in a defect-free sample. 

A solution to this problem is to display only those 
atoms whose coordination number is smaller than a cer- 
tain threshold. This number is evaluated by counting 
atoms within the nearest neighbour distance. Like the 
potential energy, in quasicrystals the coordination num- 
ber varies from atom to atom, but to a much smaller de- 
gree. As atoms near defects have a significantly lower co- 
ordination number, it becomes possible to visualise frac- 
ture surfaces and dislocation cores. For the A atoms 
the threshold for the coordination number is set to 12, 
whereas for the B atoms it is set to 14. With this method, 
the number of atoms to write out or to display can be 
reduced by three orders of magnitude. Fig. |4] shows a 
snapshot of a simulation with some 4 million atoms fil- 
tered by this technique, which was also used in a movie 
that is available online^. 

As can be seen from Fig. |4] the fracture surfaces are 
rough. Therefore to decide to which fracture surface an 
atom belongs we apply the displacement vectors between 
the initial configuration with the built-in seed crack and 
the fractured sample. The morphology of the fracture 
planes is then analysed via geometrical scanning of the 
atoms forming the surfaces. The roughness can be visu- 
alised by colour coding the height of the surface in a view 
normal to the fracture surface. 

To investigate the influence of the Bergman-type clus- 
ters on cleavage they have to be displayed together with 
the fracture surfaces intersecting them. This is done by 
restoring the initial sample without crack at zero tem- 
perature. The atoms forming the two sample parts are 
taken back to their positions in this initial sample and 



then scanned geometrically. In addition all atoms form- 
ing clusters in the vicinity of this surfaces are known. By 
displaying only these atoms and the scanned surface one 
directly can see where and to which amount clusters are 
cut. A problem of this kind of visualisation is shown in 
Fig. [5j where a cluster is cut by a flat surface. When 
looked-at from above it is obvious that only four atoms 
are separated from the rest of the cluster. When looked- 
at from below one could get the impression that the clus- 
ter is heavily intersected. On a real fracture surface clus- 
ters with centres above and below the crack surface are 
present. Therefore both views of Fig. |5] are appearing 
at the same time in a two-dimensional projection. Thus 
such pictures are not very intuitive. 

A way out of this dilemma is presented in Fig. [6j Clus- 
ters with midpoints above the crack surface are displayed 
together with the upper geometrically scanned fracture 
surface only, the other clusters are shown together with 
the lower fracture surface. As a result we get two pictures 
with clusters cut by surfaces. Note that for a qualitative 
and quantitative analysis always both pictures or sets of 
data are needed. 
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FIG. 6: Visualisation of the clusters cut by the dynamic crack. 
Midpoints of clusters are indicated as black dots, the clusters 
are idealised as spheres. The crack propagated from the left to 
the right. The upper and lower fracture surfaces are projected 
onto an X — z plane. 



5 



IV. RESULTS 

In this section the results of the numerical simulations 
are presented. For practical purposes we define k as the 
stress intensity factor K relative to the stress intensity 
factor Kc at the Griffith load: 

k = K/Kg- 

The orientations of the samples are characterised by the 
notation y,j; , where y is an axis perpendicular to the cleav- 
age plane and x is an axis in the crack propagation di- 
rection (see Fig. [6]). An axis perpendicular to a fivefold 
(5) and a twofold (2) axis is denoted pseudo-twofold (p2) 
axis. 

A. Crack velocities 
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FIG. 7: Average crack velocities for different loads and orien- 
tations. 



Simulations were performed for different orientations 
with loads in a range from fc=l.ltofc = 2 (see nota- 
tions and loads in Fig. |7]). Brittle fracture without any 
crack tip plasticity is observed irrespective of the orien- 
tation of the fracture plane. For loads below k = 1.2, 
the crack propagates only a few atomic distances, and 
then stops. Hence the energy needed for initiating crack 
propagation is about 1.4 times the value predicted by the 
Griffith criterion. Therefore, a simple global thermody- 
namic description of fracture is not applicable. The min- 
imal velocity for brittle crack propagation is about 10% 
of the shear wave velocitj^^ Vs- For loads fc > 1.2 the 
velocity increases monotonically with the applied load. 
The crack velocities are in a range of 10-45% of Vg (see 
Fig. |7|. At high loads the crack velocities on fivefold 
cleavage planes show higher average velocities than on 
the other planes. Velocities for the two different crack 
propagation directions on the fivefold planes differ sig- 
nificantly at intermediate loads (fc = 1.3). This coincides 
with ledges that are produced in the fracture surface (see 
section IV B and Fig. [s] bottom) . 



B. Fracture surfaces 

To analyse the morphology of the fracture surfaces, 
they are displayed as described in section [Til C| In Fig. [8] 
the average height is shown in grey, heights above -l-2ro 
are shown in white and heights below — 2ro are shown 
in black. The crack propagation direction is from the 
left to the right. The initial fracture surface is flat, as 
can be seen from the homogeneous regions on the left. 
The surfaces resulting from the propagation of the crack, 
however, show pronounced patterns of regions with dif- 
ferent heights. From the observation of the fracture sur- 
faces it is already evident that they are rough and that 
the peak-to-valley roughness is of the order of the di- 
ameter of the clusters. The peak-to- valley roughness 
and the root-mean-square roughness of the height pro- 
files both increas^^ for higher loads for surfaces without 



ledge formation. For fracture surfaces perpendicular to 
twofold and fivefold axes the crack fluctuates about a 
constant height (in the areas without ledges). In con- 
trast, a crack inserted perpendicular to a pseudo-twofold 
direction seems to deviate from this plane^S. 



C. Anisotropy 

As can be seen already from the fracture surfaces per- 
pendicular to a fivefold axis in Fig.[8]crack surfaces for the 
same cleavage plane differ significantly for different in- 
plane crack propagation directions. For the orientation 
5p2 ledges are produced, while no ledges form for the ori- 
entation 52 . By visual inspection of the fracture surfaces 
(see Fig. Isj) and evaluation of height-height-correlation 
functions'^ it becomes evident that for every orientation 
there are distinct angles in the height profile, which show 
pronounced height correlations that correspond to mark- 
ings in the fracture surfaces or to ledges. These angles 
are given in Table [Tj 



D. Clusters 

In Fig. |9] the clusters cut by the fracture surfaces are 
presented as described in section [III C| and Fig. |6] It is 
obvious from Fig. [9]that the dynamic crack does not per- 
fectly circumvent the clusters, but intersects them par- 
tially (right side of Fig. [9]). These intersections, however. 



TABLE I: Angles observed in the height profiles of the frac- 
ture surfaces. Angles measured clockwise from the crack prop- 
agation direction get a negative sign. 

orientation 2-2 2^ 62 5p2 p25 

angles 0°,±32° 0°,+32°,-58° 0°,±36° ±18° 0°,±90° 
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FIG. 8: Height profiles of sections of tlie simulated fracture surfaces. Load: k — 1.3, orientation 22 (top), 62 (middle), 5p2 
(bottom). The height increases from black (< — 2ro) to white (> +2ro). The scanning sphere has the same size as an atom of 
type B. 
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FIG. 9: Clusters cut by the fracture surfaces. The visualisation technique is described in Sec. |III Cj and Fig.|6] Load: k — 1.3, 
orientation: 62. 
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FIG. 10: Surface energy and density of cluster centres for the 
orientation 22. The corresponding fracture surface is shown 
in Fig. [8] top. 



FIG. 11: Surface energy and density of cluster centres for the 
orientation 62. The corresponding fracture surface is shown 
in Fig. [8] middle. 



V. DISCUSSION 



are much less frequent than for the flat seed cracks (left 
side of Fig. [9]) . More detailed analyses for different ori- 
entations validate this statement. For the orientations 
perpendicular to twofold and fivefold axes at k — 1.3 the 
ratio of clusters cut by the dynamic crack to clusters in- 
tersected by flat cuts is approximately 0.6. Additionally 
the absolute value of clusters cut by the crack for the 
fivefold surfaces is lower than for the twofold surfaces. 

Fig. [To] and Fig. [TT] display bottom up: The density of 
the cluster centres, the surface energy, a cluster in the 
proper length scale, the grey coding of the heights in 
Fig. 8] (top, middle), and the position of the seed crack 
(dashed vertical line). For the twofold fracture surface 
the low energy seed crack is located between two peaks 
in the cluster density, whereas for the fivefold surface 
this seed crack is situated close to a peak of this density. 
It is evident from the figures that it is not possible to 
circumvent all clusters by a planar cut. 

The grey coding is adjusted to the average height of the 
fracture surfaces. It is therefore evident from Fig.[TT|that 
the crack deviates for the orientation 62 from the low en- 
ergy cleavage plane of the seed crack to a parallel plane, 
reducing the number of cluster intersections dramatically 
(see also Fig. |9]and Fig. [s] middle). Samples cut flat at 
the new height show slightly higher surface energy (see 
also Fig. [nj. However, for low loads and low roughness 
the actual fracture surfaces of the dynamic cracks have 
about 5-15% higher energies than those of the low energy 
seed cracks. To assure that the dynamic crack is depart- 
ing from the initial plane not in a random manner the 
seed crack was built in at the position colour coded as 
medium grey in Fig. |11[ The resulting fracture surface 
had a similar roughness but the crack did not change to 
a parallel plane. 



Taken together the results of our simulations presented 
above indicate that the distribution of the clusters is 
crucial for the fracture behaviour: First, circumventing 
the clusters or intersecting them disturbs the propagat- 
ing crack and leads to additional radiation. This man- 
ifests itself in the crack speed. Dynamic cracks prop- 
agating in fivefold planes with few cluster intersections 
are faster than those in twofold planes, where the abso- 
lute number of cluster intersections is higher. Also the 
generation of ledges for low loads in the orientation 5p2 
costs energy and therefore slows down the crack even 
further. Second, circumvention of the clusters leads to 
characteristic height-variations. Therefore the peak-to- 
valley roughness of the fracture surfaces is of the order 
of the diameter of the clusters. Third, the observed pat- 
terns in the fracture surfaces correspond to lines along 
which the clusters are located. The associated angles are 
given by the icosahedral symmetry of the sample, namely 
0°, 18°, 31.72°, 36°, 58.28°, and 90° (see Table [!]). Ledges 
seem to be produced only for the smallest angles mea- 
sured from the crack propagation direction. Fourth, less 
clusters are intersected by the fracture surfaces than by 
the fiat seed cracks. Fifth, a seed crack at a low energy 
cleavage plane deviates to a parallel plane to reduce the 
number of cluster intersections in spite of the higher en- 
ergy required to form the fracture surfaces. In contrast, 
a crack built-in at this new position does not show such 
a deviation. 

Another observation of the simulations is that the 
plane structure of the quasicrystal also infiuences frac- 
ture. The fracture surfaces that are located perpendicu- 
lar to the twofold and fivefold symmetry axes show con- 
stant average heights. 

The three-dimensional quasicrystals give perfect cleav- 
age fracture with no indication of any dislocation activ- 



ity. This is in contrast to results on two-dimensional 
decagonal quasicrystals, where a dislocation enhanced 
fracture mechanism has been observecP^l. However a 
corresponding three-dimensional decagonal quasicrystal 
would have a periodic direction with a straight dislo- 
cation line. In the simulations presented here this di- 
rection is also quasiperiodic. As the clusters have a 
strong influence on fracture they also may bend and 
hinder dislocation lines. So dislocation emission in the 
three-dimensional icosahedral quasicrystal modelled here 
should be less likely than in the two-dimensional decago- 
nal model quasicrystal. Very high stresses are indeed 
needed to move dislocations in our model quasicrystal in 
molecular dynamics simulation^^H 

There are also indications for the stability of the 
clusters from the electronic structure of quasicrystals: 
First, experiments and ab-initio calculations show that 
directional bonding may be present within clusters of 
quasicrystals^J^^. Second, the electrons may addition- 
ally stabilise the clusters because of their hierarchical 
structur^. Therefore they should be even more stable 
than we have modelled them with simple pair potentials. 
With this evidence the results concerning the clusters 
seem reasonable and should even underestimate their sta- 
bility. 

So far fracture experiments in ultrahigh vacuum have 
only been performed on icosahedral Al-Mn-Pd, which 
has a more complicated atomic structure than the icosa- 
hedral binary model. Additionally the clusters are not 
Bergman-type. Therefore we cannot expect to represent 
this structure on an atomic level, when comparing exper- 
iments in Al-Mn-Pd with our simulations. Nevertheless, 
the size of the clusters, the icosahedral symmetry, and a 
distinct plane structure are common features and qual- 
itative aspects should be reproduced well, namely the 
size and shape of the patterns and the appearance of dis- 
tinct angles on the fracture surfaces. This is indeed the 
case, as can be seen in Fig. [T2j There a geometrically 
scanned fracture surface generated in our si mula tions is 
confronted to an STM-image of Ebert et al.'^^'^ at the 
same length scale. As we were able to correlate the ob- 
served structures to the clusters in our model, the simi- 
larities corroborate the assumption that the clusters are 
responsible for the globular structures observed in ex- 
periment. More detailed comparisons to fracture experi- 
ments in icosahedral Al-Zn-Mg-type quasicrystals would 
be desirable, but such comparative data is currently not 
available. 

We also performed fracture simulations in a C15 Laves 
phase with the same Lennard-Jones potentials to give fur- 
ther evidence, that the roughness is correlated to the clus- 
ters. This CI 5 structure is built up from deformed pro- 
late rhombohedra, one of the major tiles of the quasicrys- 
tal (see Fig. [l] and Sec. Ill A I. However, no Bergman- 
type clusters are present. When colour coded like the 
quasicrystal, fracture surfaces of the C15 cracks lack any 
roughness at low loads and cleave smoothly on high sym- 
metry plane^^. 




FIG. 12: Fracture surfaces perpendicular to twofold axes. In 
the left picture (simulation) the atomically sharp seed crack 
can be seen on the top, whereas below this area the simulated 
fracture surface appears. The orientation of the sample is 22 , 
the load was k = 1.3. The surface was geometrically scanned 
with an atom of type B. Thus atomic resolution is achieved. 
The right picture (experiment, adopted from Ebert et al.l^ 
shows an STM-image of icosahedral Al-Pd-Mn cleaved in ul- 
trahigh vacuum. As 20ro is approximately 5 nm the surfaces 
are displayed at the same length scale. 



When a crack traverses a solid, it leaves a typical non- 
equlibrium surface. Thus, up to this point, we have not 
dealt with any equilibrium surface. Even the flat Griffith 
cuts we introduced to perform the numerical experiments 
are no explicit equilibrium surfaces'^- by definition. It is 
known that i- Al-Pd-Mn surfaces sputtered and annealed 
up to about 900 K are rough with cluster-like protru- 
sions. Fivefold surf aces an nealed at higher temperatures 
exhibit flat terrace^^GSHI] which are believed to be bulk 
terminatecP^. As flvefold surfaces are not as rough as 
twofold surfaces, they are often studied in experiment. 
Adsorption on these surfaces, nevertheless, is often very 
site-sensitive. A recent review on quasicrystal surfaces 
was given by McGrath et aL^^^. All of these observations 
are consistent with our simulations of non-equilibrium 
surfaces. Fivefold surfaces experience less roughness than 
twofold surfac es (a lso in agreement with the experiments 
of Ebert et al.'^'^sl ggg pjg where more clusters are 
intersected. Additionally, since clusters are cut in the 
simulations, their binding energy is not so large as to 
term them supermolecules. Thus an nealin g at very high 
temperatures can favour flat surface^^^. Nevertheless, 
selective adsorption may then be related to completion 
of clusters. 



VI. CONCLUSIONS 

We have simulated crack propagation in an icosahe- 
dral model quasicrystal. Brittle fracture without any 
signature of dislocation emission is observed. The frac- 
ture surfaces are rough on the scale of the clusters and 
show constant average heights for orientations perpen- 
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dicular to twofold and fivefold axes. Thus both the plane 
structure and the clusters play an important role in frac- 
ture. The influence of the clusters is also seen in the 
average crack velocities for different orientations, the ob- 
served patterns in the fracture surfaces, the anisotropy 
with respect to the in-plane propagation direction, and 
the smaller amount of clusters cut by the propagating 
crack than by planar cuts. The clusters, too, are a rea- 
son why the positions of the cleavage planes cannot be 
predicted by a simple energy criterion. Since partial clus- 
ter intersections occur, the binding energy of the clusters 



is not so large as to term them supermolecules. Never- 
theless our observations clearly show that they are not 
only structural units but physical entities. 
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